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WORKING PAPER NO. 11 
THE BLANKENSHIP CONJECTURE EXAMINED 



J. D. Swift 

I August 1958 

II pages 



Hie period of x + 1 in (x? + l)/x + 1 is determined 
for $11 prunes p < 288 for which p is a prupiulve root. This 

(*$ \ 

is found to be p l 2 - 11 for all such prunes except 37* Id* 

197* 269 for which it is 1/3 of this value* The four special 
cases give counter-examples for the conjecture that the period is 
always maximal. Some arguments tending to show that the behavior is 
consistent with ’random expectation’ are given. 
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WDRKENG PAPER NO. 11 
THE BLANKENSHIP CONJECTURE EXAMINED 



J. D. Swift 

I August 1958 

II pages 

1. Introduction . For primes p having 2 as a primitive root, 
the cyclotomic polynomial 

(i) fM - 4# 

is irreducible over GF(2). It is also evident that, with respect to 
this polynomial, x has order p . The question of tne order of the 
other linear polynomial in x , y ■= x + 1 , arises. By various 
methods, one of which is included below, it is easy to see that the 
order of y is pn where n is a divisor of s = 2 m - 1 where 
m = (p - l)/2 . 

Dr. W. A. Blankinship has conjectured that n - s always. This 
conjecture was based on certain empirical evidence concerned with 
p < 100 . The chief purpose of this paper is to discuss a method by 
which the proposition was investigated for p < 288 . In particular 
the previous evidence was found to be faulty. The final results are 
that n = s for p - 3, 5, 11, 13, 19, 29, 53, 59, 61, 67, 83, 107, 
131, 139, 11 ;9, 163, 173, 179, 181, 211, 227| and that n *= s/3 for 
p = 37, 101, 197, 269. 
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Certain tables wnich were of use in the investigation and, not 
being readily available elsewhere, may be of some general interest, 
are included. 

2. Theoretical consideration s. Let the notation be as an the first 
paragraph of the introduction. Further, let g(y) = f(x) , i.e., 

(2) g(y) . £L±2 £ - L ± 

*| 

Then define z = x + “ and let h(z) *■ x" m f(x) . The degree of the 
polynomial h(z) is m . Now we maintains The order of y with 
respect to f(x) is p times the order of z with respect to h(z). 

Proof , y 2 ■ (l + x) 2 ■ 1 + x 2 » xz . The order of y 2 is the 
same as the order of y since both are certainly odd. The order of x 
is p | the order of z is prime to p . Hence the order of y is 
the order of x times the order of z by the standard theorem on the 
orders of elements on a Galois Field. Finally if h(z) divides z 11 - 1 
as a polynomial in z , m is clear that z n = 1 in the GF(2^) 
defined by f(x) . 

Thus the basic problem is reduced to the evaluation of the order of z 
with respect to h(z) or, in other terms, to finding the period 
of h(z) . 

Blankinship 1 s conjecture is equivalent to the statements h(z) 
is primitive irreducible. Now h(z) is certainly irreducible for 
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all degrees tinder consideration „ Indeed h(z) may be irreducible 
when f(x) is reducible* This is the case, for example, when p = 7 • 
The condition for reducibility of h(z) is that the corresponding 
f(x) have a proper symmetric divisor. Thus in some vague sense h(z) 
is ’more than irreducible 1 and this idea gives some credence to tne 
conjecture. It has, however, been the generally observed fact that 
there is no simple characterization of primitive polynomials any more 
than there is a simple numerical function which always yields primes. 
Indeed such functions have a statistical property known as Kronecker’s 
Hypothesis which states that the observed frequency of prunes will be 
asymptotically equal to that expected on elementary frequency considera- 
tions «, 

Now how likely is a polynomial to be primitive? The number of 
primitive polynomials is Gf(s)/m while the number of irreducible 
polynomials is 

m m 

inhere the q^ are the prime factors of m . The ratio of these 
numbers for p « 3, 5, 11, 13, 19, 29, and 37 is respectively, 1 , 1, 

1, #67, 086 , .$ho Other figures are given m a table at the end 

of the paper. Hence it is quxbe reasonable that 37 should be the 
first case of impriimtivity . Again the number of polynomials belonging 
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to e , an admissible divisor of s, is 0(e)/m . Hence if e = s/3 
the number is either 1/3 or 1/2 of the primitive polynomials -while 
if e ^ s/7 the number is ^ 1/6 of the total. Farther 3 is a 
factor of 11 of the 21 composite numbers s considered -while 7 
(whose presumed asymptotic frequency is also 1/2) is a factor of only 9 
of them. Five is never a factor. 

Hence if h(z) is impnmitive iu is most likely to have a period 
8/3 . These remarks suffice to suggest that the observed results are 

consistent with a ’purely random’ or ’Kroneeker’ behavior of h(z). 

3. The computation . To test pnmitivity it suffices to investigate 

*/<l± 

z where now the q^ are the various prJme factors of s . If one 

of these is 1 mod h(z) then h(z) is lmpramitive . Further the 

period will divide all the s/q^ which yield 1 and will not divide 

those which do not give 1 . This enables a brief calculation of the 

period. The calculation thus requires a) The polynomials h(z) for 

the required p 3 b) the prunes q. 3 c) the numbers s/q^ ; 
s/q^ 

d) sr mod h(z) . We now discuss the procedures used for these 
steps. 

a) Let f^Cx) = y?" 1 + x r *’^ + •••■»> x + 1, for r an odd 
positive integer and h^fz) » x"”^ f^Cx) for z = x + and 

k => (r - 1)/2 . Thus f and h^. are generalizations of f and n 
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to all odd and all natural numbers respectively! h m (z) « h(z) , 

The important formula iss 

(3) \(z) - s * h k-2^ * 

This recursion was first observed m a somewhat different context by 
Blankinshxp . Its proof is trivial* 

h^Cz) = x~ k (x 2k + x 2 ^ 1 + ••• + 1) by definition 

V V-1 -1 -if 

= x + x + • • • + 1 + x + * • + x 

zh. k-1 (z) + h k-2 *= (x + x“ x )x“ k+:i (x 2k “ 2 + x 2k " 3 + ••• + l) 

+ x“ k+2 (x 2k "^ + x 2 ^ + •••+!) 

«=■ x k + + (x k " 2 + + x" k+2 ) + (x k “ 2 + ••• + x“ k+2 

+ x“ k+1 + x k + (x k 2 + ••• + x" k+2 ) 

* « 

This formula gives a method of computing h(z) which is vastly 
simpler than that given by Aloert in SCAMP Working Paper 27 of 
15 February 1956. Specifically all that is neeaed is to shift 
h k _i(z) left by one and add h k _ 2 (z) . Only the output time limits 

the speed. The h k (z) 3 k < lljii were computed in less than two 
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minutes on SWAC and the specific values required Selected from the 
resulting deck* 

b) The factorization of numbers 2 11 - 1 is found in several 

tables in Kraitchiks Introduction a la Theorie aes Nombres , l-'ariB, 1952. 
Since our primes p are congruent to t 3 mod 8 (as 2 cannot be a 
quadratic residue of p) , is either odd or singly even* Hence 

the tables on pp. 12 and 38 sufficed. The factorizations are collected 
in a table appended to this paper. Hie prime factors were first placed 
on punched cards and converted to ii-preeision binary by a routine 
written for tins purpose. 

c) A division routine in Ij. -precision exact terms was written. 

This took in a number s , divided it ty a sequence of exact divisors 
and punched out the quotients. Then it accepted the next 3 . If a 
non-divisor was entered the machine halted in break -point} this feature 
guarded against typographical errors m Kraitchik or mis -punching in 
routine b). 

d) This is the principal routine and was divided into two parts. 

In the first, the input was h(z) . The routine found, by successive 
gk 

squaring, z , k = 0, 1, *•• , m , reducing the powers mod h(z) . 

As a check z =• z . The powers were stored, as produced, on succes- 
sive drum channels* The second portion accepted successively the 
numbers s/q 1 and computed z to these powers by multiplying 
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consecutively the previously computed powers which appear in the binary 
expansion of s/q^ . As a check on this routine s itself was entered 
and z = 1 computed after the maximal proper divisors had been com- 
pleted* 

The routines listed in a) s b), and c) were primarily input-output 
routines in the sense that the only time limitations were tne cyclic 
rates of these devices. The routines in d) were of rather short 
duration. The longest were for p ** 181 at 6 min. with 11 divisors 
and check and for p - 211 with 10 divisors and check, a total of 
6 min. The total run takes just over an hour for all primes less than 
288 . 

However the routine is rather hard on the machine. This seems uo 
be due to its large number of doubling commands and repeated extracts 
which cause spill and the periodic drum references following violent 
spells of computing which produce surging. It has been necessary to 
choose days of specially good machine behavior to get the routine through. 
Three such runs have been madej on these runs the single case of incon- 
sistency or failure to check occurred on s/j for p = 269 which failed 
to give 1 on the second run. 53ms particular exponent has been run 
16 times. 

As a result of these runs we can states h(z) is impmmitive for 
37, 101, 197, 269. It is highly probable that the period of h(z) for 
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these primes is s/3 • It Is highly probable that h(z) is primitive 
for all other primes p < 288 for which 2 is a primitive root. The 
difference m degrees of assertion is due to the question put* Is 
this polynomial 1 ? If the probability that the machine has run without 
error is p(m) , the probability that we should get the answer 1 by 
mistake is (1 - while the probability that we should get a 

value not 1 , shen the correct answer is 1 is (1 - p)0- - 2 -61 ) . The 
second is much greater than the first. For the three runs* the numbers 
(1 - f>Y 2~ J are so small they can be neglected entirely. The num- 
bers (l - p)^ (l - 2" m )^ s while small should not be totally forgotten. 
It must be clearly understood that in all runs mentioned the checks 
never failed j hence p is reasonably large. 

All routines are on file at NAR-TJCLA. . 
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Table 1 



P-1 

Factors of s * 2 - 1 for primes p having 

2 as a primitive root 



p Factorization of s 

3 J 1 
3 * 3 
11 * 31 
13 * 3 2 • 7 
19 s 7 • 73 
29 t 3 * h3 • 127 
37 : 3 3 • 7 • 19 * 73 
$3 i 3 • 2 731 • 8 191 
59 * 233 • 1 103 • 2 089 
61 s 3 2 • 7 • 11 • 31 • 151 • 331 

67 * 7 • 23 • 89 • 399 Ii79 

83 * 13 367 • 161; 311 333 

101 ; 3 • 11 * 31 • 231 • 601 • 1 801 • k 0£L 

107 * 6 361 • 69 1+31 ♦ 20 39b 1+01 

131 i 31 * 8 191 • lh5 293 11+3 338 ill 
139 * 7 • 1+7 * 178 l|8l • 10 032 678 938 039 
11+9 s 3 • 223 • 1 777 • 23 781 033 • 616 318 177 
163 * 7 * 73 • 2 393 • 71 H9 • 262 637 • 97 683 839 

173 t 3 * 1+31 • 9 719 • 2 099 863 • 2 932 031 007 1+03 

179 s 618 970 019 6I4.2 690 137 iil+9 362 111 

181 s 3 3 * 7 * 11 • 19 * 31 • 73 * 131 * 331 • 631 • 23 311 • 18 837 001 

197 * 3 • 1+3 • 127 • 1+ 363 933 127 297 - 1+ h32 676 798 393 

211 s 7 • 31 • 71 * 127 * 131 • 337 * 29 191 • 106 681 • 122 921 • 132 Ola 

227 s 3 391 • 23 279 • 63 993 • 1 868 369 * 1 0 66 818 132 868 207 

269 * 3 • 7 327 637 * 193 707 721 • 761 838 237 287 • 6 713 103 182 899 
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Table 2 

Polynomials h(z) for piimes of -which 2 is a primitive root 
(Notation in octal as in Marsh’s Tables of irreducible polynomials) 

3 3 

5 7 

11 67 

33 163 

19 1563 

29 73403 

37 16 3311.07 

53 7156 00067 

*59 67016 00007 

61 1 63006 00003 

67 15 63006 00003 

83 6711; 003U6 01563 

101 71 56031; 00000 33l;07 

107 670 16331; 00000 03U67 

131 67 li;030 OOOli; 00000 00003 

339 1560 31;670 00331; 00000 OOO67 

1U9 71560 OO67O 16331; 00000 03U67 

163 156 30060 00003 1;6014 00006 7H;03 

173 7H;0 33l;60 00000 06711; 0031;6 01563 

179 67340 03 US 0 00000 007H; 033l;6 00163 

181 1 633 k 0 01560 00000 00331; 07156 OOO67 

197 715 6031;0 00160 00000 00000 67OI6 00007 

211 1 56300 0711;0 33l;60 00000 00000 003l;6 01563 

227 671 1;0300 00016 30060 00000 00000 00000 711; 03 

269 7H;03 3l;600 16300 00000 03i;60 00000 00000 00000 00163 

15 Feb, 1 956, 2"^ omitted there. 



incorrect in SCAMP paper 26, 
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Table 3 

Frequencies of various classes of polynomials 



p 


m 


Irreducible 

polynomials 




Pi 


P2 


p 3 


s 




1 


— — ■ — 

i 


1 


1 




11 




6 


6 


*7S 


1 




13 




9 


6 


.S6 


.67 


.67 


19 




S6 


lt8 


.Itit 


.86 




29 


lit 


1 161 


7S6 


.283 


.65 


.93 


37 


18 


lit S3 2 


7 776 


.222 


.Sit 


.38 


S3 


26 


2 S80 79S 


1 719 900 


.ISU 


.67 


.999 


S9 


29 


18 $12 790 


18 U07 808 


.138 


.991 1 




61 


30 


3S 790 267 


17 820 000 


.133 


.So 


.33 


67 


33 


260 300 986 


231 016 608 


,121 


,81 




83 


ill 


S 3 6U7 ill SSO 


$3 630 700 7S2 


.098 


.99 96 





The third column lists the number of irreducible polynomials of degree m . 

The fourth column lists the number of primitive polynomials of degree m . 

The fifth column gives the probability that a random polynomial of 
degree m lacking a linear factor is irreducible. 

The sixth column gives the probability that a random irreducible poly- 
nomial of degree m is primitive. 

The seventh column gives (where applicable) the probability that an 
imprimitive irreducible polynomial has period 1/3 the maximum. 
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WORKING PAPER NO. 10 

THE RAND CORPORATION'S RANDOM DIGIT GENERATOR 



H. P. Edmondson 
18 August 1908 
30 pages 

The theoretical and design considerations of a machine to select 
decimal digits at random and punch them into I.B.M. bookkeeping cards 
are discussed in this report. The heart of the machine is an electronic 
binary counter which counts pulses from a random pulse source. Period- 
ically s the counter is stopped for observation. About 100,000 counts 
are expected between successive observations, so that the last digit of 
the total can be considered random. 

Analytical studies indicate that the machine is highly random in 
its selection except for trivial correlation between successful 
selections. Experimental tests of large numbers of the digits first 
tabulated by the machine Indicated no irregularities except a slight 
excess of odd over even digits. Subsequent evolution in the pulse 
forming and counting circuits appears to have entirely eliminated the 
possibility of this kind of bias. 
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WORKING PAPER NO. 15 

THE RAND CORPORATION’S RANDOM DIGIT GENERATOR 



H. P. Edmondson 
18 August 1958 
30 pages 



Introduction . Limited tables of random numbers have been published, 
but much larger tables — in fact an inexhaustible supply of random 
numbers — are needed to avoid using the same tables over and over again. 
Repetitious use of a table of random numbers is particularly undesirable 
within a single problem. 

The generation of random digit tables by human or machine methods 

1 2 

is not as simple as it appears. The remarks of Kendall and Smith * 

concerning this difficulty are pertinent: 

“It Is becoming increasingly evident that sampling left 
to the discretion of a human individual is not random, 
although he may be completely unconscious of the existence 
of bias, or indeed actively endeavoring to avoid it. 

House -to -house sampling, the sampling of crop yields, even 
ticket drawing have been found to give results widely 
divergent from expectation 

“It has long been held that mechanical methods of producing 
random series of integers do not give satisfactory results, 
pice -throwing, for example, to give a random series of the 
integers 1 to 6, notoriously results in bias. Nor are 



^M. G. Kendall and B. Babington Smith, “Randomness and Random 
Sampling Numbers, * Journal of the Royal Statistical Society, pp. l£l 
and 1^6, Vol. Cl, 1938. 

^Kenda l l and Smith, Loc. Cit., pp I52i-l56. 
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roulette tables much better, Karl Pearson has shown by 
analysis of the gaining results at Monte Carlo that the 
odds against the absence of bias are exceedingly large* 

The source of this bias is not altogether clear, but if 
we exclude the possibilities of deliberate falsifica- 
tion, it would appear to arise from small imperfections 
in the roulette wheel which direct the ball Into some 
compartments in preference to others * 

Mr. Cecil Hastings of the RAND Corporation has proposed a scheme 
for accomplishing the selection of digits with a high degree of random- 
ness, and automatically recording them at a reasonably high rate of 
speed* A machine based on a variation of his idea has been designed and 
constructed in the Development Section and put into successful operation* 
The following discussion describes Ihe operation of the machine, 
attempts to discuss its randomness analytically, and mentions a few of 
the desigi features intended to insure conformity to the theoretical 
analysis . 

Theoretical considerations. The two design criteria of a machine 

1— in I - * r , - , i 

intended to produce a table of random digits are* 

1) The device should be absolutely impartial. 

2) There should be no correlation between successive selections j 
the machine should have no memory. 

Almost any common device one might name falls down at one of these 
two criteria* A mechanical roulette wheel, for example, satisfies 
neither requirement. It is difficult to build a roulette wheel with such 
precision that one number would not be favored over another by even one 
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percent, let alone, say, one thousandth of a percent, which would be a 
more nearly acceptable figure. Furthermore, if, say, one million succes- 
sive 6*s were thrown on a roulette wheel, a groove would be worn to the 
6 compartment • Therefore, the 6 would be favored over the other 
nunfcers. 

Consider, however, the following system. Figure 1, which is a 
modified electronic roulette wheel# 




Figure 1 - Random Digit Selecting System 
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It is intended that sharp pulses from the random frequency pulse 
source should arrive at the gate at an expected frequency of about one 
hundred thousand per second* This gate circuit is controlled by broad 
constant frequency pulses from the constant frequency gating pulse 
source, so that -Hie gate allows the random pulses to pass in groups of 
about one second time duration. Each random pulse advances the position 
of the counter one digit, so that each group of random pulses advances 
the count about one hundred -thousand digits. After each group of 
pulses, the digit at -which the counter rests is considered to be random. 
This system is closely analogous to a 32 compartment roulette wheel, 
around which the ball spins about three thousand times before stopping. 

The choice of 32 numbers results, of course, from the fact that 
2^ *= 32 is the number of steps in one cycle of a five place binary 
counter. For a reason to be discussed later, the following transforma- 
tion from binary numbers to decimal digits is now used i(this trans- 
formation was not used originally). 



-U - 
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TABLE I 

TRANSFORMATION FROM BINART TO DECIMAL DIGITS 



Position in Cycle Binary Number Decimal Digit 



0 


00000 


D 


1 


60001 


1 


• 2 


00010 


2 


3 


ooon 


3 


h 


00100 


h 


5 


00101 


5 


6 


00110 


6 


7 


00111 


7 


8 


01000 


8 


9 


01001 


9 


10 


01010 


discard 


11 


01011 


discard 


12 


01100 


discard 


13 


01101 


discard 


iu 


onio 


discard 


15 


oini 


discard 


16 


10000 


discard 


17 


10001 


discard 


18 


10010 


discard 


19 


10011 


discard 


20 


10100 


discard 


21 


10101 


discard 


22 


10110 


9 


23 


10111 


8 


2 h 


11000 


7 


25 


11001 


6 


26 


noio 


5 


27 


non 


h 


28 


moo 


3 


29 


11101 


2 


30 


11110 


1 


31 


inn 


0 



The impartiality of this type machine results from the assumption 
that the pulse shaping circuit standardizes the shape of all pulses 
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It passes on to the counter, and that the level and separation of the 
standardized pulses driving the counter is sufficient to unerringly 
advance the counter one count per pulse* Even though the counter flip- 
flops themselves may prefer certain positions to others, the totals 
observed on the counter are determined entirely by -the number of pulses 
which come from the pulse forming circuits during the measured time 
intervals. It was, however, the failure of the initial circuits to 
faithfully perform these functions that caused the initial odd-even bias 
in the tables created. Impartiality also depends upon con?) let e independ- 
ence of the gating pulse generator, and the random pulse generator from 
the position of the counter. This is accomplished easily by carefully 
isolating the fields and power supplies of these different components. 

There is no evidence at present to indicate that the machine does 
not select binary numbers with complete inpartiality. However, it would 
be necessary to sample several million numbers to detect an odd-even bi$s 
of as much as one tenth percent. As insurance against the possibility 
that the machine may have an undetected bias in one of its binary 
counters, the peculiar transformation to decimal digits given an Table I 
is used. Note that the two binary numbers which transform to each decimal 
digit are complementary. Thus, if the flip-flop controlling any one binary 
place is biased by a certain amount, -die probability of any particular 
decimal digit being selected is unchanged. The excess (or shortage) in 
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the probability of the digit being selected in the first ten position^ 

1 

is exactly compensated by the shortage (or excess) in the probability of 
that digit being selected in the last ten positions. 

The effect of this complimentary combination scheme can be formulated 
analytically. Suppose that 0 ip preferred over 1 in the last binary 
place by an amount 2o( , in the next place by 2(3 , in the next place by 
28 » in the next place by 2 £ , and in the first place by 2 p . The 
probability of a 0 decimal digit equals the probability of a 00000 
binary number plus the probability of a 11111 binary number. 

p(0) - d/2 ♦ <x)(i/2 ♦ y3)(l/2 + S)(l/2 + e)(l/2 + P ) 

+ (1/2 - ot)(l/2 - pi) (l/2 - S)(l/2 - £)(l/2 - p) (1) 

Neglecting terms higher than the second degree leaves 

p(0) o l/l6 + l/8(*(3+ar£ +«& + ap + + (38 + pt + ftp + 5e +<5& + Sf + £f) 

(2) 

Similarly 

p(l) » l/l6 + l/8(-<$fi- dS -<xt -otp + /3$ + /3s + (3f + 6f + 6p) (3) 

p(2) » 1/16 + l/8(-aj8 + 96 + as + 3p - (3$ - (3£ - /Sp + St + Sf -if) (U) 
etc. 

Notice ' that no first degree errors remain as a result of this 
particular type of transformation. 
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The second criterion, the absence of correlation between successive 
selections, is certainly satisfied by this system. Actually, it would 
be nearly impossible to intentionally control the frequency of the pulse 
source and the period of the gate switching pulse closely enough that a 
"next" selection could be predicted, since the expected number of counts 
per gate interval is 100,000. The following analysis indicates how 
small this correlation actually is assuming an ideal counter, random 
pulse source, and gating system. 

The probability of exactly X random pulses occurring in any 
constant time interval group is 



P® - -g , 



(5) 



wjiere N is the expected nuoiber of pulses per group. 

If, therefore, the count starts from a digit d Q , the probability 
of its advancing just k digits to digit is 



pOO 



ic r 



-N 



( 6 ) 



The digit d^ Would also be selected if the counter advanced 
32 + k counts, and the probability of this happening is 

^32 ^ 

" (k + 32} 1 
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p(k + 32) 



(7) 
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Similarly, the digit can be selected by the count advancing k 
plus any multiple of 32 counts. Thus, the entire probability of the 
digit d^ being selected after d Q is 

„ fJ , _ N 3 * _-N . N^ 32 _-N N** 61 * _-N . ... , 0 \ 

H e + T FVl ' gyf e + (FV15C)T ® * * (8) 

Simplification of this to a finite series can be achieved by the 
use of the identity 



Ijk ^32 j,k *6h 

ki * wT'wn * trmryi + * 4 



hi 



± n 

31 ikg- Ne 



e 



e 



(9) 



m=o 



Thus, 



p(V - ^ 




31 

£ 



m=o 






(10) 



This equation reduces easily to the form. 

( , i -n S- Ncos tf“ l(N TS"‘ k fr) 

P\°k) * 52 e E e 6 

m=o 



( 11 ) 



Since N is about 100,000 , the term in the summation corresponding 

i 

to m = 0 is by far the most important. Next are the two terms corres- 
ponding to m «• 1 and m « 31 , and the remaining terms are negligible 
in comparison with these. The three retained terms can be written 
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1 -N 



pCdfe) - e 



N Nc0{ $ i(Nsin^-k^) 



e + e 



+ e 



"“■T 1 ( SslD T5 L - k Tr> 

q e 



• (12) 



But cos jg * cos > sin ^ «* - sin J and. for k an integer, 

ir* 31 if 

k radians is coincident with -k radians . Thus 

if" -N(l-cosg)| i(Nsin^-kJk) -i(NAin^“ k^)\ 1 

pCdjj.) <= -jg L 1 + e l e +e 

-N(l- cos^j) 



1 

3$ 



1 + 2 e 



cos (N sm ^ - k y&) 



(13) 



Since the cosine function can be no greater in absolute magnitude than 
unity, then this probability can differ from the perfect value of 1/32 
by no more than 



|l/32 - ^ e 



-n(i ■ 



cos ^) 



(Hi) 



For N ■ 100,000 this deviation is 



|l/32 -p(d k )l ^e- 100 ^ 1 - 980 ^ 




(15) 
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Indeed, the correlation between successive selections is negligible. 

The random pulse source . The circuitry used as a random pulse 
source is a high gain wide band noise amplifier followed ty a detector 
biased so that only the noise peaks above a certain high level are trans- 
mitted through the detector into the output circuit. Figure 2 is a 
schematic of the circuits used. The source of random noise is simply 
shot effect in the first vacuum tube. The r-m-s value of this noise is 
controlled by the bias applied at the input grid. The overall bandwidth 
of the amplifier is about 6 megacycles. 

The justification for using a highly biased random noise detector 
as a random pulse source may not meet the approval of the critical 
reader. However, all that is needed from a practical standpoint is a 
highly irregular and unpredictable source of pulses to drive the counter, 
and the biased random noise detector certainly satisfies this requirement. 

As a matter of fact, it can be argued that the pulses generated by 
such a device are nearly truly random. The requirement of a truly 
random source would be that the probability of a pulse occurring between 
t and t + dt should be some p dt , where p is the expected 
number of pulses per second and is a constant entirely independent of 
the number and distribution of pulses generated up to time t . 

Figure 3 is a typical random noise voltage signal, with one 
detected pulse shown to aid in discussing the problem. 
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Suppose that at time t Q a pulse is finally observed. Then, however, 
the observer is able to predict a trend for a short interval ahead. Know- 
ing the intrinsic decay behavior of the amplifier an question, he knows 
that this decay voltage superimposed on thp new random signal voltage 
increases (or decreases as the case may be) the probability of a detect- 
able peak being observed. After At, however, the decay trend will 
have expended itself, and the probability of a pulse will remain constant 
(so far as the observer knows) until another pulse is observed. 

The length of At can be assumed to be less than one microsecond 

for the amplifier in question, since one-twelfth microsecond is the 

conventional rule-of-thumb decay time estimate for a low-pass amplifier 

of six megacycles bandwidth. The fact that the amplifier is actually 

band-pass instead of low-pass can be neglected, since the ratio of noise 

power in the missing low end of the frequency range to the power in the 

band pass region is quite small. It will become apparent later that this 

* 

divergence from pure random occurrence in an interval of one microsecond 
following each observed pulse is of no extra concern* The pulse forming 
circuits reject any pulse which falls within one microsecond of a pre- 
viously observed pulse anyhow. 

The gate circuit and the gating pulse generator . The function of 

r 

the gate circuit and its controlling gating pulse generator is to 
measure out intervals of one second during which pulses from the random 
pulse generator are amplified and passed on to the pulse shaping circuits. 

- lit - 
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Between each of these intervals should be a period of about one -tenth 
second during which the pulses from the random pulse generator are 
blocked, and the digit at which the counter stops is read and recorded 
in an I.B.M. card. This timing sequence is obtained easily ty means of 
an unsymmetrical multivibrator. 

Figure is a schematic of this timing multivibrator and the gate 
circuit. The gate circuit is ^ two stage pulse amplifier with the plate 
of the second amplifier tube tied in common with the plate of the gating 
pulse isolation tube. Notice that when the multivibrator lies with its 
right-hand tube conducting that the isolation tube is cut off. Thus, 
the pulse amplifier works as a simple amplifier with no interference from 
the multivibrator isolation tube. When the right-hand tube is cut off, 
however, this isolation tube grid goes positive with respect to its -90 
volt cathode, pulling its plate down to a negative value. Consequently, 
Hie voltage is removed from the plate of the second pulse amplifier 
tube, and no pulses from the random pulse source can pass to trip the 
pulse forming circuits which follow. 

A second multivibrator isolation tube is shown dn Figure it. When 
the grid of this tube goes positive, its plate current actuates the 
counter reading relays and, subsequently, the I.B.M. key punch. 

The pulse shaping circuits . The pulse shaping circuits have a 
difficult job to perform. The Input pulses are of various sizes and 
shapes and occur at random in time. From these highly irregular input 
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signals, the pulse shaping 1 circuits must form output pulses of standard 
size and shape, no two of which can be within one microsecond of each 
other. Any input pulse which cranes within one microsecond of a previous 
pulse must be rejected. 

Although this would ordinarily be an easy function, it is made very 
difficult by the rigid requirements placed upon the dependability of the 
circuits. The pulse shaping circuits must be perfect, lest counter 
partiality bias the digit selecting system. If any pulse leaks through 
the pu2,se shaping circuits which is of such a size and shape, or is so 
dose (less than one microsecond for the counter used) to another pulse, 
that the counter might or might not (according to its own preference) 
advance one count, then the unavoidable partiality of the counter itself 
is allowed to contribute partiality to the system as a whole. The 
circuits of Figure E> were arrived at experimentally, and appear to he 
absolutely dependable. 
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First in the circuit are two cascaded one-shot multivibrator 
circuits (6SN7 , s), with a pulse width of about one microsecond* These 
circuits do the major portion of the work. Pulses which are too weak to 
trip the circuits do not get through to the output at all, and nearly 
eveiy pulse strong enough to trip the circuit produces a standard one 
microsecond pulse. Some complications, however, occur. For example, when 
two strong pulses occur just about one microsecond apart, the second 
pulse may catch the one-shot multivibrator circuit just as it is resetting 
and produce something different from the standard pulse shape. 

The final insurance against irregularities is a relatively fast flip- 
flop circuit of 6AK5 ’s driven at the cathode of one of the 6ak£’s by 
a 6L6 cathode follower. This circuit trips to the right if the driven 
cathode is made more positive than +10 volts, and resets to the left if 
the driven cathode is made more negative than -10 volts. Thus, the 
circuit can trip only once for each full pulse from the one-shot multi- 
vibrators, and any small input irregularities will fail to trip the 6AKJ> 
flip-flop unless they go from -10 to +10 volts. There is little 
chance of trouble with input pulses too closely spaced m this circuit, 
because the flip-flop is veiy fast compared to the circuits which drive 
it. The flip-flop has a rise time at its plates of about one tenth micro- 
second, and stabilizes in an exchanged position at the end of about two 
tenths of a microsecond. 
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Thus, the output of this circuit is a square wave of about one 
microsecond duration and very steep sides. If this output is differen- 
tiated (not shown in Figure 5 ee figure 6) the result will be 
gharp positive and negative pulses, and, obviously, no two positive pulses 
can be closer together than one microSepond. Also, these pulses will be 
of" standard-' site ahd" Shape, s±nb6. the relatively sluggish driving circuits 
ahead of the 6AKS flip-flop cannot effect its rise time appreciably# 

The counter circuit . The function of the counter circuit is to 
accurately count the pulses as they come from the pulse shaping circuits* 
The resolving time of the counter must be short compared to the one 
microsecond minimum spacing between successive input pulses, and the 
circuits must be absolutely dependable else counter partiality might 
contra. bite partiality to the system. 

The counting is done in the binary number system because it is the 
natural system for electronic counters. All the binary places except 
the last five are disregarded, giving a count cycle of 2-* *= 32 steps, 
the 33rd step being identical with the 1st, the 3l*th step being 
identical with the 2nd, etc., etc. 

In Figure 6, the first two tubes are a cathode follower to isolate 
the flip-flop of Figure 5, and a driver tube for the following flip-flop. 
Recall that the output of the flip-flop in Figure 5 is a square wave 
pulse with very steep sides. 
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In the output of the cathode follower circuit this voltage is differen- 
tiated by the 5>0 pji condenser into the 10,000 ohm resistor, so that 
the input to the 6AK6 driver tube is a sharp positive pulse correspond- 
ing to the front edge of the input pulse, and a sharp negative pulse 
corresponding to the trailing edge of the input pulse* Of 'course, the 
positive pulse only is effective, since the tube is normally biased 
beyond cutoff. When this sharp positive pulse is applied to the driver 
tube grid, its plate conducts momentarily. 

The basis of the electronic counter is a d-c flip-flop, a circuit 
having two stable positions. Note m Figure 6 that if no external 
driving pulses are supplied, the 6AK£ flip-flop would sit with either 
its right hand tube conducting and the left hand tube cut off, or with 
its left hand tube conducting and the right hand tube cut off. Note also 
that if the grid of the 6AK6 driving tube (which is normally cut off) 
is pulsed with a sharp positive peak causing its plate to conduct 
momentarily that both plates of the 6AK5 flip-flop will be momentarily 
brought to almost zero voltage, and the flip-flop circuit will then 
return to the state opposite the one it was resting in when the driving 
pulse occurred. This exchange of position is caused by the two ■memory" 
condensers shown (20 jjp). The condenser to the formerly non-conducting 
plate has a greater voltage across it than the condenser to the conduct- 
ing plate* Thus, if both plates are momentarily reduced to nearly zero 
volts, the condenser having the greater voltage across it causes the 
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grid of the opposite tube to be the more negative. When released, then, 
the opposite tube becomes non-conducting, while the formerly non- 
conducting tube conducts. Each positive pulse from the pulse shaping 
circuits reverses the position of the first flip-flop. 

The output of this flip-flop, then, is alternate positive and 
negative steps. This output is isolated by a cathode follower and 
differentiated by circuits similar to the differentiating circuit used 
ahead of the first flip-flop. The result, of course, is again vezy 
sharp positive and negative pulses, and the positive pulses are used to 
drive the second flip-flop. This chain — flip-flop, cathode follower 
stage and driving tube — is repeated a total of five times. The first 
flip-flop reverses conduction tubes for each input pulse. The second 
flip-flop in the chain reverses for each positive pulse from the differ- 
entiating circuit following the first flip-flop — and this occurs only 
on evezy second input pulse. Similarly, the third flip-flop reverses 
on evezy fourth input pulse, etc., etc., and the last flip-flop reverses 
on every sixteenth pulse. Thus, the counter cycle consists of the 
following 32 steps* 
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The binary to decimal transformation . Table I of the section on 
theoretical considerations shows the manner in which the positions of 
the binary cycle are to be interpreted as decimal digits. It is 
desirable to indicate each of the digits by a closed circuit rather than 
by a light, voltage, or current, so tha^ any type of automatic device 
such as an electric typewriter or I.B.M. duplicating punch may be used 
to record the selections. Figure 7 shows how this may be accomplished 
ly the use of multi -pole double-throw relays. 

Note that, depending upon what combination an which the relays are 
open or closed, any one — but only one — path is closed to the common 
input point. If these relays are controlled ty -the position of the d-c 
flip-flops of the binary counter chain, then it is possible to determine 
by inspection Hie particular combination which closes the circuit to each 
particular output point. In Figure 7 each of the 32 output points is 
labeled with this binaiy counter combination (assuming 1 to mean upper 
contacts and 0 to mean lower contacts), and the decimal digit this 
combination should represent is copied from Table I* Then the two out- 
put points that indicate each of the decimal digits are tied to a 
common output terminal. 

figure 7 shows how the relays a ip driven by thyratrons, the gmds 
of which are controlled by the cathode follower voltage of each of the 
five flip-flops. The gating pulse output tube m the lower right hand 
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comer of Figure it closes the master relay - power relay in Figure 7, 
famishing plate voltage to the five transformation relays. Thus, the 
transformation relays do not attempt to follow the progress of the 
electronic counter, hut merely are controlled by the counter during the 
one -tenth second interval during which the gate is shut and no pulses 
are driving the counter. 

The relay -power relays work in conjunction with the master 
transformation circuit relay to prevent the application of power to the 
recording circuit until all five transformation relays have been given 
ample time to set in the selected combination, and to open the recording 
circuit before the transformation relays are released. Otherwise, each 
time the gating pulse occurs, the transformation relays would give 
momentary false circuits when they were pulling in or releasing. 
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Test and monitoring circuits . Unfortunately, it is impossible to 
devise any kind of monitoring circuit which will indicate whether or not 
the machine is choosipg numbers without bias. There are, however, a few 
simple circuits which chedf the performance of the components of the 
system. 

Referring to Figure U, note the "Test Sequence Switch" in the grid 
circuit of the 6AC7 gated pulse amplifier. When this switch is closed, 
normal bias is applied to the grid of the 6AC7, and the circuit func- 
tions as a pulse amplifier. If the switch is open, however, a bias of 
-10£> volts is applied to the 6AC7 grid, and the tube is completely cut 
off* Thus, no pulses from the random pulse generator can drive the 
counter. However, every tune the gating pulse clamps this 6AC7 circuit, 
one pulse is jFbrmed by the pulse shaping circuits pf Figure J? (the lead- 
ing edge of th,e gating pulse being sharp enough to trip the pulse 
shaping multivibrators). ULth this happening, the counter should advance 
just one count per selection. This test sequence is valuable as a 
check on the reliability of the system froifi the gate circuit through 
to the output device (which is an I.B.M. key punch at the present time). 
Before and after each running peripd the machine is set on "test 
sequence" for several minutes, and the punched cards produced are 
checikpd for errors. 

A second circuit monitors the average 1 rate at which pulses are 
being counted. This circuit is a simple electronic frequency meter 
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connected to the last flip-flop in the counter chain. A voltage is 
generated proportional to the average frequency of the counts and is 
indicated by a meter on the front panel of the instrument (see Figure 8). 
Also, this voltage is used as an a-v-c source to regulate the average 
random pulse rate to the frequency desired (about 100,000 per second) 
by using it as a "Control Bias" on the noise source tube in Fig. 2* 

Electro -mechanicjal counters were installed in the output circuits 
to indicate total counts. Ten counters were connected directly across 
the output terminals of the transformation relays to indicate the total 
number of times each decimal digit is selected. Also, two counters were 
installed to count the number of times the first flip-flop of the 
electronic binary counter indicated right and left (0 or 1) — a measure 
of the impartiality of the system up to the binary selections. Recall 
that the binary-to -decimal digit transformation used. Table 'I, yields 
decimal digits of improved impartiality. Thus, it is advisable to look 
for partiality in the binary selections since partiality would be more 
evident there. 

Conclusions ♦ 

1« A machine which takes as random the last digit of the total 
random pulses in a fixed period has been constructed and pub into 
successful operation. Initially, an improbable excess of even over odd 
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selections occurred. But, since revisions were made in the original 
circuits, no indication of partiality has occurred. 

2. A complimentary type of transformation from binary to decimal 
digits is used. If partiality should exist m the binary numbers 
selected, this particular type of transformation would yield decimal 
digits with considerably less partiality. 

3* The machine is completely automatic. The unit built was 
connected to an I.B.H. key punch ,to compile a table of several million 
random digits. 
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WORKING PAPER No. 20 

A CLASS OF MAPPINGS AND SOME EXAMPLES 

G. A. Hedlund 
28 August 1958 

This note concerns a class of mappings defined and 
studied by Rothaus (SCAMP Working Paper No. 25, August 30, 1957) 
and Blackwell (Ibid and SCAMP Working Paper No. 3, July 9, 1957) . 
The first part of the paper develops some of the general theory. 
Most of the theorems proved are not new, though it appears that 
the proofs are. The second part of the paper gives some examples 
which seem to have been unknown previously. 

Let S^ denote the set of all sequences of 0's and I's of 

length n, n a positive integer. Any member of S n will be called 
an n-block. 

Let S denote the set of all unending sequences of 
O’s and l’s . Any member of S is a function with domain 
the set I of all integers (positive, negative or zero) and range 
m Sj, . Let 

s - .... s^s^s^... 

and 

t * .... t ^ f 0 ^*1 c • * 

be members of S. We define a distance in S as follows* 

d(s,s) = 0 . 
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If t |s , there exists a least non-negative integer k such 
that s^ £ t^ ors^^t^., and we define 

d<s - t) * dr ■ 

It is easily verified that S , with this metric, is a Cantor 
discontmuum. 

Let f be a function with domain S n and range m 
Sj , i.e. , f * S n ~> Sj . Then f defines a mapping g m of 

S , mto S ,as follows; Let B c S . and let 
m+n-1 m m+n-i 

B = s lSz s m+n -l . 

Define 

t^ = f(s^, • • ® * S x+ n _l^ * l = coo, m, 

and let C a t, t 9 t . Then g (B) * C» 

Similarly, f defines a mapping g of S into S, as follows 

Let s c S and let s = . „ . s , s s, s, ,,, . 

-I O l Z 

Define 

tj S f (S i* S i+1»*’** S i+n-l* ’ 1 e 1 ' 

and let 

t — • « • t ^ t^ ■ , a # 

We define g (s) - t. Clearly g is continuous, 

A basic problem is to determine conditions on the 
function of f which will assure that the mapping g is an 
onto mapping, 1 , e, , 



g (S) = S 




Remark. 



g is an onto mapping if jand only if, 

g m ^m+n-1^ * for all m , or, equivalently, if, and only if, 

B being any m-block, there exists an (m+n-l)-block C such 

that g (C) * B . 
m 

The stated condition is obviously necessary. 

To prove the sufficiency, suppose that, B being any 

m-block, there exists an (m+n-l)-block C such that g (C) - B 0 

m 

But then g (S) is clearly dense m S . Since S is compact 
and g is continuous, g (S) is compact, thus closed, and, being 



dense m S, g (S) = S. 

Le t s s . . o s i s s,... 

-1 o 1 



belong to S. Then s is said 



to be periodic if there exists a positive integer p such that 

s . = s , 1 e I . 

l+p l * 

The least such positive integer is the period of s . 



Remark. 



If s « S is periodic, then g (s) is periodic. 



If s has period U then the period of g (s) divides a) „ 

In the following if A is a set, crd (A) denotes the 

number of members of A . If B c S , g~* (B) denotes 

m “m 

the collection of all members C of S , such that 

m+n-i 

g m (C) * B. In general, if D is a subset of S m , g^ 



(D) denotes the collection of all members E of S 



m+n-1 



such that g (C) e D , 
m 
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Lemma. 



Let the mapping g defined by f be an onto 



mapping and let there exist a positive integer k such that 
crd g m ^ (B) * k for some m-block B and crd gp^ (A)z k 

for all p-blocks A and all positive integers p. Then 
crd gp^^ (BA) = k for all p-blocks A and all positive 



integers p. 
Proof. 



It is sufficient to prove that 

crd < B0 > * k * crd *m« (B1) 



Let 

*m < B > - tC l , C 2 C k ] . 

A member of g m +j (BO) cannot be identical with a member 
of g m “J (Bl) . Thus crd g m “j [BO , B1 ] 21 2k . 



But 

g m+l t B0 * B1 l C [ C i°- C 2° C^O.Cjl.^i,...,^ 1] 

# 

and thus 

g m+l [BO.Bl] = [CjO.CgO C^.C^.C^ Cy] 

and crd g m+ j [BO, Bl] = 2k . 

Since crd g m ^[ (BO) z k, crd g m ^J (Bl) z k , 



it follows that 

crd g m+l < B0 > * k * crd g m+l * B1 >* 

The conclusion of the lemma now follows by induction. 



- 4 - 




REF ID : A38876 



Theorem 



Let the mapping g be an onto mapping. Then 



— 1 n— 1 

crd g (A) ^ 2 for all p-blocks A and all positive 
integers p. 

Proof. Let I* denote the set of positive integers and let 

k a Min crd g * (B ) 

m e r\ B c S 
* m m 

Since g is an onto mapping, kZl . Let B bean m-block 
such that crd g | (B) a k. It follows from the preceding 
lemma that crd g (BA) a k for all p-blocks A and all 



positive integers p. 
Let 

-1 



gm (B) * t°l' C 2 C k ] * 



Then 



*m+q (B V C t C l S q * .... C k S q ], qcl + . 

We recall that S denotes the set of all q- blocks and BS 
q M q 

denotes the set of all (m+q)-blocks with initial block B. Let 

the collection |C,S , . C, S be donoted by C 

1 1 q* * k q J 7 q 

♦ q 

The set C has k * 2 members. 

q 

Now different members of BS cannot have the same 

q 

inverses uqder g m ^ q and each member of BS q has exactly 

k inverses under g . ^ . Thus crd g (BS ) ** k« 2^ , 

m+q m+q ' q 



«m+q (BS q ) S C q and the mapping 8 m+q <C*) a BS q is 



exactly k to 1 



- 5 - 




REF ID:A3887 



Let A be an arbitrary p-block and let crd g p (A) = w. 

Let q = p+n-1 . Now S contains exactly w members, the 

$ 

image of each of which under g is A„ Thus m C there are 

& p q 

exactly kw members whose images under g m ^.q are blocks ending 

1 

in A. Now BS contains exactly 2 ~ blocks ending m A. Since 

* 

the mapping g m +q* BS is exactly k to 1, the image set 

* 

under g , of kw members of C must contain at least w 
e m+q q 

n—1 

members. It follows that 2 ” >. w and the proof is completed. 



Theorem. 



Let the mapping g be an onto mapping. Then 



- 1 n- 1 

crd g p (A) a Z ~ for all p-blocks A and all positive integers p. 



Proof. 



Let A be a p-block and suppose 
crd g" 1 (A) £ 2 n ~ 1 



- 1 n - 1 

From the preceding theorem we infer that crd g" (A) < 2 and 

"1 n—1 ^ 

crd gq (B)< 2 for all q-blocks B and all positive integers 

q . Now g (S .) * S , crd S , , “ 2P +n- * and 

^ & p p+n-1 p p+n-1 

crd S p * 2 p . Also crd [S p+n-1 - g p X (A)] > 2 P+n " 1 - 2 n ' l (2 P -1) , 



g [S . , - g" 1 (A)l bS - A 

e p 1 p+n-1 6 p ' p 



e p 1 p+n-1 e p ' p * 

1 1 

Not more than 2 members of S . . - g” (A) can 

p+n-1 °p 

map, under g p , into the same element of S p - A . 



N(B, qi)< 2 qi - 2 qi • 2" n+1 
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Thus 



2 P - 1 = crd [S p - A] > 2" n+1 crd [S p+n _ 1 - g' 1 (A)] > 2 P - 1. 

From this contradiction, we infer the truth of the theorem. 

Remark, The preceding theorem shows that the property 

that g be an onto mapping is equivalent to the property that 
g be noisy m the sense that it transforms a random sequence 
(all blocks equidistnbuted) into a random sequence. 

Theorem, Let g be an onto mapping and let s « S . Then 
crd g 1 (s) < 2 n_1 . 

Proof. Suppose that crd g * (s) > 2 n ~* . Let 



-1 



a. I 

g (s) as tj, t^, . . o , t^. . Consider any pair t^t , i f j . 

There exists an integer p^ such that the central (2p^ + l)-blocks 

of t and t are not identical. But then the central (2p+l)-blocks 
^ J 

of t ,t differ for all p > p . Let p be an integer such that 

i J y 



p > max fp | 1 < i < j <k] 

Then no two of the central (2p+l) -blocks of , t 2 



t • • • » 



*k 



are alike. But the images under g^ p these k blocks 

are identical. If k > 2 n * , this contradicts a preceding theorem. 

The proof of the theorem is completed. 

Remark. It appears that the multiplicities of the mapping g 

at different points may be different. It would be of interest to 
investigate the various possibilities and characterize them. 
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If g is an onto mapping and 

S “ m c SiS S-|S 0 «ee 

-1 O 1 C 

is a periodic sequence, then each member of g ^ (s) is 

periodic. Let s have period to , let t c g"^ (s) and let p be 

n— 1 

the period of t . Then p * p to , and 1 < p < 2 
Proof. Let g be an onto mapping, let 

s * • • * ® 

-1 o 1 2 

be periodic with period to and let t e g ^ (s) . Let p (k) be 

s 

the number of different k-blocks in s . Then p g (k) < to for 

all k By a preceding theorem we infer that t contains at 
11“ 1 

most 2 63 different (k+n-l)-blocks for each positive 

integer k. Thus, if p^. (m) denotes the number of different 

n _l 

m-blocks m t , we have P t (m) < 2 03 for all m . 

Suppose t has no period less than 2 + 1 From 
lemma 7. 2 (Morse and Hedlund^ American Journal of Mathematics, 
Vol. 60, 1938, pp 815-866) p^. (m) > m + 1 for all values of 

m for which p^ (m) < 2 n + 1 . But this is true for all 

m and thus p^ (m) >. m + 1 for all m Let m « 2 n-1 63 
1 1 

Then p fc (2 to ) > 2 60 + 1. contradictory to 

IX* 1 

P t (m) < 2 to for all m . We infer that t has a period 
less than 2 n +1 and t is periodic. 
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Let p. be the period of t . Then p. < 2 n ”^ w From 

a preceding theorem, there exists a positive integer p such 

n- 1 n— l 

that p. * p co , and we have o> < pco < 2 and 1 < p < 

Lemma Let m be a positive integer and let B be an 

m-block. For q _> m , let N(B, q) be the number of q- blocks 
which contain B . Then 






q °° 



Proof. 



We first observe that N(B,q)/2^ is a montonic 



increasing function of q. For if C 1 s a q-block which contains B, 
then CO and Cl are different (q+1)- blocks each of which 
contains B . Thus 



and consequently 



N(B , q+1) > 2 N(B, q) 

N(B , q+1) ^ N(B,q) 

2 q + l ^ -pj— 



Thus we can assume that q = pm and it is sufficient to 
prove that 



lim 

P 00 



N(B , pm) 
2 pm 



= 1 



Let B = B^ , and let B^, B 2 B^, k = 2 m , 

be the set of all m-blocks. Any block of length, pm can be 
written as a p-block of B^s. There are k^ such blocks. Of 
these there are (k-lp which do not contain B g B^, and 
thus k^ 5 - (k-1)^ which do contain B. Thus 
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and hence 



N(B , pm) > k P - (k-lf * 2 mp - (2™ -if 



and 



N(B,pm) > i - n . 1 p 
_pm — ,m r 



But 



lim 

p — ^ oo 



(1 ~±f . 0 

2” 



lim N(B,pm) 



p — > oo 



jPm 



= 1 



The proof is completed. 

Lemma . Let B be a block. For q > n , let D denote 

a partition of the set of all q-blocks into sets of q-blocks each 

■n— 1 

containing Z members. For q sufficiently large, all members 
of some element of the partition must contain B as a 
sub- block. 

Proof. We suppose the theorem false. That is, there 

exists a sequence of integers q, < q_ < . . . , D ,D and 

1 * q l q 2 

partitions such that some member of each element of fails to 

contain B as a sub-block. But then the number of members of 

D which do not contain B is at least equal to the number of 
qi 

elements of D , or 2 qi /2 n- * . Thus, using the notation of 

the preceding lemma, we have 
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N(B,qi) 






< 1 



rr 



This contradicts the preceding lemma. 
Definition* 1 The sequence 

s * 



®-l S o S 1 8 2 



is said to be transitive provided every finite block appears in s . 
Theorem. Let g be an onto mapping and let s be transitive. 

Then each member of g *(s ) is transitive. 

Proof. Let B be an arbitrary k-block. The collection 

of all (m+n-l)-blocks 



-1 



g (A) I A e S m defines a partition D 



m+n-1 



,n-l 



into sets each containing 2 blocks . From the preceding lemma 

we infer that for m sufficiently large there exists an m-block 

A such that each member of g ^ (A) contains B . Now A 

°m 

appears m s and each member of g ^ (s) must contain an 



-1 



element of g (A) 



-1 



It follows that each member of g (s) 
contains B and thus is transitive. 

Let f be a function with domain S and range in S, 

and let , g^ n ^ be corresponding mappings of S m+n ^ 

into S and S into S , respectively. Let f be a function 
m p 



with 



domain and range m , and let g m 
the corresponding mappings. 



(P) „(P> 



g 



be 
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The mappings g|^+ p _i and can be composed in 

an obvious fashion to define mappings g ;** 1 • g' of 

m m+p-1 

S , , - into S and thus a mapping • g^ of S into S . 

m+n+p-2 m 

Similarly there is defined a mapping g^ » g^ of S into S . 

It is not necessarily true that g^ * g^ * * g^ . 

Let s t S 

8 s ooe 8 1 S 0 S L © ° ° 



and let 



t * o o o t ^ t^ tj e o o 

be defined by 

*i * s l+ i • 1 * 1 • 

The transformation <j> s s — ^ t is called the shift transformation. 
It is a homeomorphism of S onto S whose properties have 
been studied extensively (see Gottschalk and Hedlund, Topological 
Dynamics, Am. Math. Soc . Colloquium Publications , vol„ 36, 1955, 
Gh. 12). A subset Y of S is invariant if <j> (Y) - Y . 

It is easily shown that the transformation g of S into S, 
defined by f , commutes with i.e. f g <{> s <|> g o 
Theorem. Let f define an onto mapping g of S onto S and 



and let x be a proper closed invariant subset of S , Then 
g (X) is a proper closed invariant subset of S. 

Proof . Suppose g is an onto mapping and X isa 

proper closed invariant subset of S. Then X is compact, 
g (X) is compact and g (X) is closed. Since g $ * <f> g, 
g (X) is invariant. 
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Suppose g (X) s S. Let s be a transitive point and 
let x « X such that g (x) = s . According to a preceding 
theorem, x must be transitive and thus X * S , contrary to 
hypothesis. The theorem is proved. 

Corollary. g^*t g^ n ^ is an onto mapping if and only if both 
g ^ and g are onto mappings. 



Proof. 



Clearly, if g^ and g^ are both onto mappings, 
the g' ' is an onto mapping. 

Suppose g^ g (n *(S) C g^(S) £ contrary to the 
supposition. Thus, in any case, g^* ^ is onto. Now if g^ 

In) 

is not onto, then g' '(S) is a proper closed invariant subset. 
Hence, by the last theorem g^ g^(S) is a proper closed 
subset of S , again contradicting the assumption that 

1S onto. Thus g^ must also be onto and the second 
part of the corollary is proved. 

The remainder of this paper is devoted to the determination 
of all functions f which determine onto mappings g' ' for 
n/4. 

The totality of functions f which define onto mappings 
m the cases n s 2 or 3 are easily compiled and are as 
follows. For n s 2 there are 6 such functions of which 
three are 



f{x L , x 2 , ) 



= x. 



= X. 



■ *1 + x 2 
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and the other three are the duals of these, that is, the 



functions obtained by adding 


1 to each of the function values 


For n * 3 there are 


30 such functions of which 15 


as follows: 


f(x lt x 2 ,x 3 ) 


1 


X 1 


2 


X 2 


3 


x 3 


4 


x l + x 2 


5 


Xj + x 3 


6 


*2 + *3 


7 


Xj + x 2 + x 3 


8 


Xj + x 2 x 3 


9 


x 3 + x l x 2 


10 


xi + x 2 + x 2 x 3 


11 


x i + x 3 + x i x 2 


12 


x L + x 3 + x 2 x 3 


13 


x 2 + x 3 + x L x 2 


14 


x i + x 2 + x 3 + x L x 2 


15 


x l + x 2 + x 3 + X 2 X 3 



and the other 15 are the duals of these. 

Remark. Of the fifteen listed, the first six are compositions 

of mappings for which n s 2 „ 
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For the case n = 4 , it is considerably more difficult 
to determine which functions determine onto mappings. It is 
known that if f is linear in either x^ or , that is, f 



is defined by 



f(x 1 ,x 2 ,x 3 ,x 4 ) =x x + f 1 (x 2 ,x 3 ,x 4 ) 



f(x 1 ,x 2 ,x 3 ,x 4 ) =x 4 + f 4 (x 1 ,x 2 ,x 3 ) , 
then the corresponding mapping is onto. There are 496 such 
functions 

It is also known that if f is defined by composing a 
pair of mappings of lower order (m this case a 3 and a 2) then 
f defines an onto mapping if and only if each of the composing 
mappings is onto. It is easily verified that there are 22 such 
composed functions which are not linear in x^ or x 4 and 
which defme onto mappings Of these 11 are given m the following 



table 



f(x 1 ,x 2 ,x 3 ,x 4 ) 



x 2 + x 3 
x 2 + x 3 x 4 

x 3 + X 1 X 2 

xj + x 3 + x l x 2 

x 2 + x 3 + x l x 2 

x 2 + x 3 + X 3 X 4 

x, + x. + x_x . 

2 4 3 4 

*1 + x 2 + X 3 + X 1 X 2 

x 2 f x 3 + x 4 + x 3 x 4 
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and the remaining 11 are the duals of these. 

Now a necessary and sufficient condition that the mapping 

g defined by the function f be an onto mapping is that 
-1 n-l 

crd (B) p 2 for each m-block B, and each positive 

integer m . But it is sufficient (theorem due to Blackwell. 

see Rothaus, loc.cit. ) that crcl g~* (B) a 2 n ~* for m * 2 n 

m 

n- 1 

and each 2 -block B. 



This criterion is not as difficult to apply as first 



appears if use is made of the following device, illustrated for 



the case n s 4. 

Let the set of all possible 3-blocks be denoted as follows* 



000 0 

001 1 

010 2 

(1) 011 3 

100 4 

101 5 

110 6 

111 7 



Then the 4-blocks can be denoted 



0000 


00 


0001 


01 


0010 


1 2 


0011 


1 3 


0100 


24 


0101 


25 


0110 


36 


0111 


37 


1000 


40 


1001 


41 


1010 


52 


1011 


53 


1100 


64 


1101 


65 


1110 


76 


1111 


77 
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XX ( 1 , j*0, 1, 2, .... 7, j s 2i or 2i + 1, mod 8) 

i J 

denotes the 4-block of which the initial 3 -block is x and the 

l 

terminal 3-block is x 

J 

Now for any specified function f (on the 4-blocks) we 



list under 


0 , those blocks B for which f (B) * 0, 


under 1, 


the complementary set. 


For example: 




_0_ 


_ 1 _ 




00 


12 




01 


24 




13 


37 


(3) 


25 


41 




36 


52 




40 


64 




53 


65 




76 


77 



Any 5-block B can be written in the form abc , where 
a, b and c are integers from 0 to 7 , a is the integer corresponding 
by (1) to the initial 3-block of B , b is the integer corresponding 
by (1) to the middle 3-block of B and c is the integer corresponding 
by (1) to the terminal 3-block of B . Thus 01011 can be written 
as 253 „ 

Now if the 5-block B a abc is to map (under g^) into 00 , 
then ab and be must appear under 0 in (3) and this condition 
is clearly sufficient. It is thus possible to obtain the 5-blocks 
which map into 00 by taking any element ab under 0 in (3) 
for which the second term b appears as a first term and following 
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ab by the second term of all elements for which b is the first 
term. As an illustration we determine the 5-blocks which map 
into 00 and 01 respectively 



00 


01 


000 


012 


001 


137 


013 


252 


136 


364 


253 


365 


400 


537 


401 


764 


536 


765 



The process can be continued in an obvious fashion. The 
7 -blocks which map into 0000 are 

0000 

00000 

00001 

00013 

00136 

40000 

40001 
4001 3 
40136 

But it should be noted that in continuing this process, it is only 
the digits (representing 3-blocks) which appear at the ends which 
are of concern to us and the intermediate ones can be suppressed. 

We list, in terms of their terminal 3-blocks only, the blocks which 
map into 0, 00, 0000, 00000000 , respectively. 
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_0_ 


00 


0000 


00000 


00 


00 


00 


00 


01 


01 


01 


01 


13 


03 


03 


03 


25 


16 


06 


06 


36 


23 


40 


40 


40 


40 


41 


41 


53 


41 


43 


43 


76 


56 


46 


46 



When it was found (by hand computation) that there 
existed a function (n»4) which defined an onto mapping, which 
was not linear in any variable and which was not obtained by 
composition of onto mappings of lower order, it appeared that 
it might be worthwhile to determine all such functions. 

This determination was carried through on SWAC. The 
non- restrictive assumption was made that f(0, 0,0, 0) « 0 . 
Then of all such functions, those were rejected which did not 
produce an equal number of 0's and l's (eight of each). The 
remainder were then successively subjected to tests as to 
whether crd (B = 8 for each 2-blocks , crd 
g~^ (B^) a 8 for each 4-block , crd gg* (Bg) = 8 for 



each 8-block Bg . There were 291 such successful matriculants. 
The following data concerning these was recorded on punch cards. 
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1) A tabulation of the function f . 

2) The eight 4-blocks constituting f (0) m terms 
of their terminal 3-blocks. 

3) The coefficients of the polynomial defining f . 
The machine time in carrying through this program 

was 80 minutes. 

Of the 291 functions defining onto maps, the 248 linear 
in the first or last variable were sorted out, leaving a residue 
of 43. Of these 11 were known to be obtainable by compositions 
of lower order onto mappings, leaving a residue of 32 „ These 
3 2 functions are tabulated m the following two tables; 
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Table I Functions on 4-blocks which determine onto mappings, which are not linear in 
either the first or last variables and which are not obtainable by compositions of lower 
order (n < 3) onto mappings 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 

0000 00000000000000000000000000000000 

0001 000000000000000000000000111111 11 

0010 00001111111111111111111100110101 

0011 00110011111111111111111101010110 

0100 11111100001111110011101011111111 

0101 11110100110001110101100100110000 

0110 11010011010010011000011111000011 

0111 11001111100100001110010010101001 

1000 00000000000000001111111100000000 

1001 010110110010011001100011000011 1 1 J, 

1010 10101011111110000001100000111011 N 

1011 001100001111111110000100010101 10' 

1100 11111101101010110011010111111111 

1101 10101101101010110011010111000000 

1110 01110110010101001100101011001100 

1111 11001110010101001100101010101001 
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„ % 






\ 



Table II. Non-zero coefficients of the polynomials defining the 
functions listed in Table I. 

f(Xi* X2* 5(3* X4) " ^ a i ^1 ^ Xj Xj ^ C ijk ^1 ^k 

l<i<4 1< i<j<4 l<i<j<k<4 



*2 



12 



B 



13 



J 14 B 23 ®24 



34 



C 123 C 124 C 134 ^34 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 1 

18 1 

19 1 

20 1 
21 1 
22 1 

23 1 

24 1 

25 

26 

27 

28 

29 

30 

31 

32 



1 

1 

1 

1 

1 1 
1 1 
1 
1 
1 
1 

1 1 
1 1 
1 1 
1 1 
1 1 
1 1 
1 
1 

1 1 
1 1 
1 1 
1 

1 1 
1 

1 

1 

1 1 
1 1 
1 

1 1 
1 

1 1 



1 

1 

1 

1 

1 

1 

1 

1 



1 

1 

1 

1 

1 




1 

1 

1 

1 




1 

1 

1 1 
1 



1 1 
1 
1 
1 

1 1 
1 
1 



1 1 
1 

1 

1 

1 1 
1 

1 

1 

1 1 
1 
1 




1 

1 

1 

1 

1 

1 

1 

1 



1 

1 

1 

1 

1 

1 

1 

1 



1 

1 



1 



1 
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It was pointed oat by R A Dean and R. C Lyndon 
that the set o£ 32 polynomials listed in Table II, can be 
generated from a set of eight by application of two simple 
processes One of these is the substitution (x^x^) 1 e„, 

interchange of x^ and x^ -and interchange x^ and x^ 

The other is complementation, 1 e„, substitution of 1 + x^ for 

x.^ , 1 =1, 2, 3, 4 The following is a generating set 



7 


x 3 + *4 


10 


x 3 + x 2 x 4 


24 


Xj^ + X 3 + X. 


22o 


Xj + X 3 + X. 


9 


x 3 + x 1 x 2 


13. 


x 2 + x 3 + x 


14 


x^ + x 3 + X 


23 


x^ + x^ + X 
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WORKING PAPER No. 22 

A NUMERICAL SIMULATION TO STUDY THE READINGS OF A 

FIXED ANTENNA ADCOCK RADIO DIRECTION FINDER 

UNDER CONDITIONS OF FADING 

Joseph F. Mount 
C. Tompkins 
25 Feb 1958 
10 pages 

In this paper we shall describe a numerical simulation 
from which we try to draw some conclusions concerning the 
behavior of a perfectly balanced fixed antenna Adcock direction 
findei [1, 2, 3] under conditions of fading. It is assumed that 
fading is caused by random cancellation of signals arriving 
on different transmission paths. Normally these paths might 
differ either in the number of reflections from the ionosphere or 
in the layer from which they are reflected, or both. Minor 
variations in directions of the transmission path from any 
jump pattern may occur at the receiving site because of 
variations in ionospheric tilt, variations in penetration, local 
reflections, and other causes. 




to 
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UN CLASSIFIED 

will be centered about two sets o£ direction cosines, which 
may be specified independently of each other and of other 
simulations at the start of each calculation; an approximate 
intensity is specified precisely for each transmission, inde- 
pendently of all other inputs. We allow for the following random 
variations in transmissions, each variation to occur slowly 
enough to justify treatment through use of variables which 
are piecewise constant: 

(a) The relative phase between the two arrivals 
varies randomly and uniformly between 

0 and 2ir ; 

(b) The intensity of each arrival will vary from the 
approximate value set by a factor 1 + h, where 
h is a sample from an approximately normal 
distribution with mean zero and standard deviation 
which is set as part of the input of the problem; 

(c) The six direction cosines of the arriving signals 
will be modified to direction numbers by the 
addition of independent samples from normal 
distributions each with mean zero and with 
standard deviations which are set as part of . 
the input to the calculation. 
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These eight standard deviations are set independently of 
each other and independently of other parts of the calculations. 

At least for the time being the calculation will consist 
of the determination of two functions and the recording of 
relations between them. For each choice of the random variates 
the signal strength received at one of the collectors will be 
computed and for each choice the reading of a fixed antenna 
Adcock direction finder will be estimated. The signal strength 
will be calculated straightforwardly. The reading of the 
radio direction finder will be estimated either by computing the 
points of maximum deflection of a Watson Watt [2] type display 
under the stimulus of the two signals, or by computing angles 
of minimum deflection for a rotating goniometer [l] type 
display. In each case, it will be assumed that the equipment 
receiving the signals is m perfect adjustment, so that none 
of the calibrating complexities of [l] and [2] will be introduced 
except to the extent that the angle of arrival influences the 
octantal error. 

It is proposed that relations between these two functions 
be exhibited by means of tallies m one hundred positions, 
giving a kind of scatter diagram. Thus, if the functions are 
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denoted by f and F , we propose that the range of possible 
values of each be divided into ten parts, arbitrarily as part 
of the input routine, and that the number of examples found 
m which the functions fall into each possible pair of cells be 
recorded and presented as output. 

2. The functions to be computed -- the two-channel case. 

For the two-channel case, we examine equations (8) 
of [l] (these equations will be repeated below), and try to 
compute the response of a two-channel direction finder 
under the influence of a pair of signals of the type indicated 
above . 

Explicitly, we assume that signals with frequency 
u)/Zis are being received over two paths. The quantities 
pertinent to the two signals will be distinguished by subscripts, 
the subscript 1 for the first signal and the subscript 2 for 
the second. The first signal will arrive with direction cosines 
a ^ , pj and , the second will have direction cosines 

a Z ’ ^2 ' ant * ^2 ° The S1 8 na ^ s wl ^ ^ve amplitudes (weighted 
by the effectiveness of the collectors) and A^ respectively. 

The first signal will be taken to set the standard in phase, 
so that the responses of the two antenna pairs to it will be 
given by direct analogues of equations (8) in [l] , 
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( 1 ) 



wPjd 

E NS1 = 2A 1 Sin ~Zc~ 



Gj<*jd 

e ewi * **1 aln -2T- 



s m cot , 



sin cot , 



where, as before, d is the distance between collectors 
of the pair and c is the velocity of radio wave propogation. 
For the second arrival the phase is displaced a constant 
amount sampled from a uniform distribution over the 

range [0, Z v ] . (Here plays a role somewhat different 

from the role of £ in [ 2 j. ) The responses to the second 
signal are respectively 



E 



<op 2 d 

NS 2 = ^2 Sm ~ZT Sin (wt ‘ 



(2) 



COOT 2 d 

E EW2 = ^2 Sm -2 — Sin <“ l - % Z ) 



The analysis here follows that of [2] 



The total NS signal is the algebraic sum of the two 
signals received? 

w{3jd o>{3 d 

E N g - 2[A^ sin — 2 — sin wt + A sm — g sin (wt- f^)] , 
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or 

u&d up_d 

(3) E ns * 2[(A J sin— 2 ^- + A £ cos % z sin — — ) sin cot 



w P 2 d 

- A^ sin ^ sln — 2c — COS 01,4 ■* “ 

Similarly, 

(oarjd u>(3 2 d 

(4) e ew S 2 L( A i sin— 2 ^ — + A 2 cos § 2 sin — ^ — ) sin 



ut 



A_ sin §■ 



sin 



“P 2 d 

2c 



cos ut] . 



The development of the functions to be tallied follows 
the general development in [2] 0 

A heterodyning signal cos (u + u Q )t is mixed 

with E^g and E^,^ , the difference frequency is extracted, 

and the signals amplified to be presented respectively as 
horizontal and vertical deflections. Using (3) and (4) and 
the relations 

cos (u + u )t sm cot * ^r[sin (2w + u )t - sm u t] , 

(5) ° 

cos (u + U )t COS ut a -tfcos (2u + u )t + cos u t] 
the retained signals are proportional to 
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CO 



'Ey = -(A.J sm- 



f*i d . . . "P 2 d 



+ A cos fj_ sm 



~Tc~ ' “2 *2 2c 



) sin u t 
o 



(6) 



<op 2 d 

- A_ sint sin = — 

2 ?c 



cos to t 
o 



(jor^d ojo«_a 

\E tt = -(A, sin — = + A_ cos £, sin — = ) sm u> t 

H '1 2c 2 2 2c o 



r 2 = E v 2 + E 2 



wp d 

A, sm sm — = cos w t 

2 2 2 c o 

The deflection may now be calculated as a function of t* 
formally. 

The maximum deflection (and the minimum) occurs when 
dr/dt « 0, hence when 

(7) E y E V + E h = 0 * 

where the dot denotes differentiation with respect to time. 

Since we propose a computational, and not an analytic, 

attack on die problem, we propose to simplify the notation 

at this point. To that end write the computable functions 

(op.d wp_d 

1 « A > 



V„ a 



S a - (Aj sin-jj- + A 2 cos e 2 sm — ^ 

wPzd 



) . 



( 8 ) 



a- A 2 sine 2 sm 



uar,d 

Hg a - (A^ sin + cos ^ 



uott-d 

2 sm -2r-> 



H c a - A 2 sin § 2 sm 



*<* 2 d 

2c 
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so that 



(9) 



Ey 


s v s 


sm 


(0 t + 
o 


V C 


cos 


(0 t 
o 


e h 


s H s 


sin 


(0 t + 
o 


H c 


cos 


(O t . 

o 



Then the condition (7) for maximum or minimum 
deflection gives 
{10 ) tan 2 



. . „ 2 (v c v s + H c H s» 

“°° V c 2 - v| + H 2 - H s 2 



In a complete cycle of 2 ir radians, this equation (10) is 
satisfied by four angles (o Q t o separated by ir/2 radians. 
Values for Ey and E^ a t these points may be 



found by writing 



-1 ± Jte 



tan co t s 
o o 



tan 2w t +1 
o o 

tan 2(o t 
o o 



using a similar formula to find tan ( 0 Q t o /2 (choosing the 
ambiguous sign carefully each time) and writing 



sm to t 

o o 



( 11 ) 



2 tan (o t /2 

o o' 

_ - 

1 + tan (o t /2 
o o' 



1 - tan (o t / 2 
. o o' 

COS (O t SB , - 



o o 



1 + (o t f 2. 

o o' 



Independently of the choice of sign, the four dis- 
placements may be computed easily. If the computed values 
of (11) are written 
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S = sm til t and C = cos w t , 
o o o o o o 

then at the four extreme deflections the components are 

E,/ = S V Q + C V_ and E 1 > S H_ + C H , 
V 0S0C H oo o C 



( 12 ) 



'E ,. 2 = C V c - S V_ and E 2 * C H Q - S H_ . 

V 0S0C H oS oC 

E 3 =--S V 0 - C V,, and E 3 = -S H - C H-. , 

V oSoC H oS°C 



,E V 4 * - C o V S + S o V C and E H - - C c» H S + S o H C • 



3 4 12 

Ey and Ey above are obtainable from Ey and Ey by 



multiplying by -1 , and a similar remark applies to the 
horizontal components, only the first and second of these 
deflections will be considered. Write 

(E®) 2 *= (Ey®) 2 + (E H ®) 2 , a = 1, 2 , 



the first or the second of the deflections will be retained 

12 2 2 

depending on whether (E ) or (E ) is larger. We 
choose the first function to be tallied as 

F L * max ff [(Ey®) 2 + (E H ®) 2 ] . 

is a measure of the received signal strength. 
Finally note that for the deflection retained the dis- 
placement is at an angle <t> o where 
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(13) 



E 



tan c)> 



Of 

H 



a 



For some a , |3 write tan (j) Q = 



for the chosen displacement index a . This is the tangent 
of the angle read on the direction finding equipment 

ya . Choose 

a , p to be horizontal direction cosines accepted as correct-- 
say the mean from which and p^ deviate. Then the 
tangent of the "error "in reading (including octantal error) is 



tan <j> - tan cj> 

* 14 * F 2 = tan ^ " V *1 + tan <|> tan $ ‘ 

o 

This function F^ is the second function to be tallied. 
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